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I. INTRODUCTION 


Considerable interest has recently been shown in trap- 
ped waves travelling along the boundaries of continents. A 
"waveguide" effect exists over the continental shelf. That 
is, wave energy is confined (essentially by refraction) to 
the continental shelf. Two general types of these waves 
exist: 

n» Edgewaves which are characterized by wavelengths of hun- 
dreds of kilometers (km) and periods of hours (almost always 
less than a pendulum day). 

B. Shelf waves, which are generally even longer, have pe- 
riods greater than one pendulum day, and travel southward 
along the west coast of an ocean in the northern hemisphere 
(as do Kelvin waves). 

STOKES (1846) showed that such a wave guide effect ex- 
ists over a uniformly sloping beach or continental shelf, 
with the amplitude of the gravity waves decaying exponen- 
tially to seaward. URSELL (1952) showed that Stoke's edge- 
waves were the fundamental mode of a family of waves order- 
ed by the number of modes parallel to the coast. REID(1958) 
studied long waves on uniformly sloping shelves of infinite 
width, including the effect of the Coriolis force. Reid 
showed that the sea surface may react as an "inverse baro- 
meter" and that atmospheric pressure systems may be a driv- 
ing force for edgewaves. He found that the Coriolis force 


could cause the wave period to vary from 46% less than to 


ES 


86% greater than that for the non-rotating case, depending 
on the direction of travel. A new quasigeostrophic wave is 
now possible, analogous to a Kelvin wave, having.no small 
scale counterpart. 

ROBINSON (1964)initiateda study of the continental 
shelf wave and studied the data of HAMON (1962, 1963) re- 
lating tidal and barometric conditions at several stations 
on the eastern and western coasts of Australia. In this 
model the continental shelf ends abruptly, at which point 
the depth becomes infinite. He found that an inverse baro- 
meter effect was exhibited but that the propagated shelf 
waves had a celerity double that of his calculations for 
the western boundary. MOOERS and SMITH (1968) studied the 
relation of sea level and barometric conditions along the 
Oregon coast for a period of nearly one year. Their statis- 
tical results show a barometric factor of -1.2 cm/mb and 
predominant sea level oscillations of 0.1 and 0.35 cycles 
per day in the summer. They conclude that a shelf wave of 
period three days is travelling north. MYSAK and HAMON 
(1969) found shelf waves off the coast of North Carolina, 
but found no coupling between the sea surface and air pres- 
sure in the frequency range 0-0.5.cpd. ADAMS and BUCHWALD 
(1969) show that an equally suitable driving force for 
shelf waves is the longshore component of the geostrophic 
wind. This may account for the exaggerated frequency re- 
sponse of the sea level on the east coast observed by 


Hamon. 
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MYSAK (1967, 1968) extended Robinson's work, and discus- 
sed the effect of a continental shelf of finite width on the 
frequency of Hamon's Australian waves. His theoretical so- 
lutions correspond more closely to the observations, al- 
though he still cannot account for the extremely low read- 
ings along the eastern boundary. He attributes the discrep- 
ancy mainly to the presence of stratified water and currents 
in the deep water beyond the continental shelf. A signifi- 
cant discrepancy exists between the dispersion relation and 
that for waves over an infinitely wide continental shelf. 

This paper is a study of the effects of a continental 
slope and finite ocean depth upon the present one-slope 
models of MYSAK (1968). A sharp discontinuity in the depth 
of water beyond the continental shelf is not a common occur- 
rence in the world ocean. It is interesting to study the 
two-slope situation where a gently sloping continental shelf 
(slope, s < 0.002) and steeper continental slope (s e 0.05) 
form a transition zone between the coast-line and deep 
water. Three parameters: the slope of the continental 
shelf, the slope of the continental slope, and the depth of 
the deep water should have possible effects on shelf waves. 


These are investigated below. 
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II, ANALYSIS 


= = 
h s x (x W.) 


1 
9 RA 
Й tÀ = = 
“چ‎ \ i i h, SX d Wı (Ss; 52) 
! h | 
y na (W, x = молат) 
ha 
ga j | h,= constant 
=9 W, mn, де Е 
ІІ 1 ΠΠ sS 8ے + ۷ے‎ 


Е19 а. Ргоб ето Еце"Ссма 7s Lope smoda lg 


The model is characterized by a gradually sloping con- 
tinental shelf of finite width (Region I) adjoining a steep- 
er continental slope (Region II) which terminates in water 
of uniform depth (Region III). A representative slope and 
width of the continental shelf of .002 and 100 km (Mysak, 
1968a) are used below. A representative depth of the 
deep ocean is 5000m and is the greatest depth of Region III. 
Three slopes will be used for the continental slope: .03, 
.05 and .08, with .05 used as the standard for comparison 
with Mysak's model. 


The shallow water equations are used: 
oc tye π 
du/ot - fv + g “fox = dvfot + fu+g °/ду=0 (1) 


o/dx(hu) + d/dy(hv) + d€/dt = O (2) 
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where (u, v) are the (x, y) velocity components and ё 


is the free surface height. 


Consider a wave, moving in the y direction, speci- 


fied by 
u, = u, (x) ei (987m) 
vw, 7 V, Go) el (007m) (3) 
" E 


where k denotes the region. 


Using (1) and (3), u and v are 








u = I (fmn - σ -- j ei(ct-my) (4) 
f^-g k 
αμα een i(ct-my) 
MO MI, یت‎ ME (5) 


Eliminating u, v from (2) in Region I gives the equa- 
tion 


28 -2 s. fm 





T i o -f l _ 2 = 
hymno. s m SO sS KS a یراول نو‎ (6) 
After making the substitutions 
hı = SX 
_ ا‎ 7 im 
Р] gs; σ 
(6) can be written 
xT, + т + (p, mə>x) n, = 6 (7) 


را 





This equation in 7 has as its solution (REID, 1958) 


1 


-2,/2 
пу = 1A,Fa, (2,) + B,Ga,(z,) Je (8) 


where Fa, (2,) is Kummer's function, and Ga, (2,) is a sec- 
ond (independent) solution (SLATER, 1960): 


2 3 
Fa (z) = 1 + аг + Sete , alatiplaraje +... 


n 
R a mem ο. (9) 


(n:) 





and 
1 ET δε ή 
Galz)  ے‎ 7 ۷ ہت‎ Т 
nal 1 0 2 
+ а(а+1) ••• (а+п-1) 2 A + η E 
(10) 
m-p; 
Here, а, nn Zi — 2mx, and Άη» В, are constants 


of integration. Since Ga4 (zi) approaches infinity as Z 


approaches zero, Ва must equal zero: 


-2,/2 | 
n = A, Fa, (z,)e -(11) 


Substituting (11) into (4) gives 


i E dFa; (zi) 
Ui = = 2 а A, q (£+0) Fa, (2,)-20 uw Ж (12) 


Similarly, ingpegion ΤΙ, 


-2,/2 
п = е (Аға. (2.,) т Bea (13) 
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and 








-тЕ/? 
; 2 dGa.(z.) - 
_ _ Ägme_ { [ В _ Be? 
U, dme. A, |(0-£) Fa,(z,)- 20 — ] 
O -f 2 
= аса 
+ B_| (0-£)Ga_(z_)- 20— 252! 7 (14) 
اد‎ 57 dz, 
where 
Mae 
25 = 2055, а, = ,کس‎ 
h 2 
m A 2 
2 Ξ апре gs, с? 


In Region III the counterpart of (6) is 


DM? 
и c^ -f 2 . 
h,n, + ( = - ham ) n, = S (15) 


Since N must be bounded for large x, 


па = деви where (16) 
: 2 "m M 
fne | m^ - وا ےہ8‎ and 
да 
3 
igA, 2 (x-W, -W,) 
U. = - —— )070+0 (٥ه‎ (17) 
3 2 


There are now equations defining and U in each re- 
gion. The next step is to patch together the solutions for 


n and U at the points x = W ande (W, +W.) thus eliminat- 


1% 
ing the constants Aye By, . The patching conditions are 
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2 B 
LE = E = 0 (surface height continuity) (18) 
2 E 
Chul, = [πα], = (normal flux continuity) (19) 
where 
1 
ا‎ eS 


The following abbreviations will be used: 


ds = 0-2.7 + = za) 

The subscript, 1, refers to the solution for Region I 
(continental shelf) where it joins Region II (continental 
slope). The subscript, 2, refers to the Region II where Te 
joins Region I. The subscript, 3, refers to the continental 
slope where it joins Region III, the flat bottom. The func- 
tions  subscripted 3 have the same form as those subscripted 
2 with the exception of the variable, 25, which is deter- 
mined by the distance from the origin. 

Using (18) and (19) between Regions I and II and set- 
ting A. = l, the constants of integration A, and B, can be 


| 


solved for: 


G.F.' = F.G.' 

tt (20)‏ ك = رھ 
G,F, - F4G,‏ 
,ھ2 PF.‏ 

в, = il 2 2 L us 
Got 2 
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Using (18) and (19) between Regions II and III gives the 


final equation in terms of m and © only: 


GNE - РС , 
¿A yotm) E, T 2m, ' | 
+0 πα 
Ὃς 2 22 
It υπ | 
ка 21 (ше, + пе, - 
ο Εώς Е СО 


(22) 


Because there is no way to solve (22) analytically, it 


is necessary to find the roots numerically using the IBM 


360/67 computer system at the Naval Postgraduate School. 


For a fixed w = O/f a search routine was used to find the 


several m's satisfying the equation. The computer work is 


described in Appendix A. 
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ТТТ. RESULTS AND CONCLUSIONS 


A. REVIEW 

There are a number of questions to be answered. MYSAK 
(1968a) in his finite-width model found: 

1. Shelf-wave numbers are inversely related to the 
shelf width, for a fixed frequency. 

2. There is a low wave number cut-off for edgewaves 
which is a function of the shelf width. That is, as the 
shelf width increases, the smallest possible wave number 
decreases (the largest possible wavelength increases). 

3. The fundamental mode edgewaves of REID (1958) 
with periods greater than one pendulum day do not exist 
over a continental shelf of finite width. 

It should be noted that Mysak's solution is only ap- 
proximate in that the maximum shelf depth h is assumed much 
less than the ocean depth, D, leading to an approximation 
of the equation expressing continuity at the edge of the 
shelf (Mysak's equation (10)). His results (labeled MA 
below) depend on this approximation. Exact solutions of 
Mysak's equation (6) (corresponding to equation (22) in this 
work), were also generated so that comparisons could be 
made among MA, an exact solution for Mysak!s model (ME) and 


results of the two-slope model studied in Section II(TS). 
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m. COMPARISON OF MYSAK'S APPROXTMATE SOLUTION WITH AN 
EXACT SOLUTION FOR MYSAK'S MODEL 
Nine cases were studied in order to compare ME and MA. 
Two of the cases are illustrated in Fig. 2. In all cases, 
the continental shelf is 100 km wide with a slope of .002, 
duplicating Mysak's sample calculation. The bottom depth 


varied from 500m to 5000m. The shelf wave results are 


shown in Figs. 5 - 9. 
Lata! πα "E ^ Ru 
DLE | πηι". 1 

ы | a Z 

В. --. 4 y | 

i Л а 

7 (000 rr E 3c00 rn 

|-—/00 km — w 7777777 O س‎ , E 3% 


(Vertical exaggeration = 100:1) 


Fig. 2. Profile of the finite-width model. 


Note that: 

1. The approximate solution consistently gives a 
smaller wave number for a particular w and mode. It ap- 
pears to be the limiting condition for the exact solution. 

2. Except for the fundamental mode, an error of less 
than 1% exists between corresponding modes of MA and ME in 
cases where the depth ratio h/D is smaller than .067. 

3. For the fundamental mode there is still a signifi- 
cant error introduced in m when using MA for large depth 


ratios and frequencies (>10% for w> 0.8 when D=5000m). 
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4. The edgewave results are shown in Figs. 18 - 20. 
Neither MA nor ME gives a fundamental edgewave mode similar 
to that of Reid. This can be seen directly from Mysak's 
equation (10) in the approximate case. Because the argu- 
ment is positive by definition, a must be negative in order 
for the Laguerre function to have roots (zeroes). This in 
turn requires that either f > G or я >> f. 

5. A low wave number cut-off does exist for ME edge- 
waves, which diminishes with an increase in depth and in- 
creases with an increase of | ω | (Table I). The cut-offs 
are always lower than those for MA, and are quite symmetric 


with respect to direction of travel. 


Table I. ME edgewave wave number cut-off (mW) 


Deep-water depth (m) First mode Second mode 
1000 008 1.32 
3000 Оса 0.76 
5000 0726 0.58 


As pointed out by Mysak and Reid, the edgewave disper- 
sion relation is not symmetrical with respect to direction 


of travel, due to the influence of the Coriolis parameter. 
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©, COMPARISON OF THE MYSAK MODEL WITH THE TWO-SLOPE MODEL 
۱ /00Am | ka < /G0£m —=| 56Ка F«— 
lokm | رت‎ .᾿ M. 









2 | 
DA κ 3000 со 


d A 


(Vertical exaggeration = 100:1) 


Pee. 3. Two-slope model. 


The Mysak model was compared with a two-slope model 
whose continental slope was .05 (Fig. 3). The results are 
shown in Figs. 5 - 9. Note that: 

1. Except for the fundamental mode for small & ,there 
is little similarity between the dispersion relations for 
the two models, especially for the large deep-water depths 
(i.e., wide continental slopes). For a deep-water depth of 
5000 m, both mode 2 and 3 waves of TS have smaller wave 
numbers than mode 2 of ME for any given frequency. 

2. The fundamental TS shelf wave does not asympto- 
tically approach w= 1.0 for large wave numbers as does the 
corresponding ME wave. In fact, the fundamental wave now 
behaves like the fundamental edgewave of Reid. 

3. A low wave number cutoff is still present for 


edgewaves, but at a significantly lower wave number (for 
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mode l, nearly half that of ME). The cutoffs are no longer 


symmetric with respect to direction of travel (Table II). 


Table II Two-slope edgewave wave number cutoff (mw, ) 





Deep-water depth ш < 0 w> 0 
πι Mode 1 Mode 2 Mode 1 Mode 2 
1000 0.42 132 0.28 123 
3000 0.26 0.76 0.16 0. 46 
5000 009 0.58 012 058 


Do COMPARISON OF CONTINENTAL SLOPES 





100 km — w | - A Ge 77 
333 Ke. 1212 km 


Slope = 0,03 Slope = 0.08 
(Vertical exaggeration = 100:1) 
Fig. 4.  Two-slope models with different continental 


slopes. 


Three values were used for the value of the continen- 
tal slope: .03, .05, and .08 (Fig. 4). The results are 
shown in Figs. 10 - 18, for deep-water depths from 500 m 


to 5000 m. 
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1. Wave numbers for a particular mode of trapped 
waves over a fixed deep-water depth decrease with an in- 
crease of continental slope width. 

2. For the fundamental mode over a constant gradient 
continental slope, wave numbers increase with an increase 
in slope width (i.e., an increase in deep-water depth). 
All other modes decrease with an increase in width. 

3. A curve is not available for mode 3 for a conti- 
nental slope of 0.03 and deep-water depth of 5000 m. This 
15 attributed to unknown problems of the computer routine. 
This problem does not occur elsewhere. 

4. Discontinuities appear in the dispersion relations 
for modes 2 and 3 with a continental slope 0.05 (in the 
deep-water depth range 2200-3400 m) and with a continental 
Slope 0.08 (for depths greater than 2800 m), suggesting the 
presence of complex values of ш. A similar phenomenon is 
not observed for a slope 0.03. Equation (22) was investi- 
gated for complex values of m and real w (Appendix B) and 
00 roots were found for a slope of 0.05 and depth 
2800 m (Fig. 22). Since the surface height is the real 
part of Č (x, y, t) complex values of m imply a spatial 
growth rate exp (Im (my) ] in the positive y direction. The 
roots m are complex conjugates, so that one wave grows and 
one decays at this rate. Then the most unstable wave 
(i.e., the one with the maximum growth rate) would be ex- 


pected to dominate the shelf-wave spectrum. 
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5. No complex wave numbers were found for the edge- 
waves studied. Continental slope width is inversely pro- 


portional to wave number as in Mysak's results. 


E RECOMMENDED FURTHER STUDY 

The next step would be to study further the effect of 
different continental shelf slopes using the TS "c й 
More study is mandatory in D 4 above, both to: 

l. rind its ба ES e 

2. determine if this is a mathematical curiosity or 
a physical reality. 

Future investigators should seek to avoid approxima- 
tions to their models. As this study has shown, the results 
for the exact equations can be significantly different than 


the approximations. 
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APPENDIX A 
NUMERICAL SOLUTION OF ONE-AND TWO-SLOPE MODELS 
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